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We unveil the novel physical origin of the insulating block checkerboard antiferromagnetism in 
vacancy-ordered K0.8Fe1.6Se2. Our first-principles electronic structure analysis reveals its incom- 
patibility with a simple Fermi-surface nesting or Mott insulator scenario, and suggests the picture 
of coexisting itinerant and localized electronic states. Consistently, we demonstrate that it can 
be unified with the metallic collinear or bicollinear antiferromagnetism of the vacancy-free parent 
compounds LaOFeAs, BaFe2As2, or FeTe in the spin-fermion model. These results indicate that 
the blocking effects of Hund's rule coupling and the resulting electron correlation are crucial to the 
electronic and magnetic structures of iron-based superconductors. 

PACS numbers: 74.70.Xa,75.10.-b,71.27.+a,75.25.Dk 



One of the puzzling phenomena in iron-based super- 
conductors (FeSCs) is that despite apparent similarity 
in crystal and electronic structures, their parent com- 
pounds exhibit different metallic antiferromagnetic (AF) 
patterns: collinear C-type (Fig. la) and bicollinear E- 
type (Fig. lb) in pnictides [11 and chalcogenides [3], 
respectively. Moreover, a transition between the C and 
E types in the same material can be induced by merely 
varying the anion height from the Fe plane [4^ . Such mag- 
netic softness implies the presence of strong spin fluctu- 
ations and electronic correlation, which are generally be- 
lieved to be at the heart of the high-Tc mechanism [5j |6] . 
It is thus urgent to resolve this puzzle and to classify the 
essential nature of electronic correlation in FeSCs. 

Recently, a new horizon to look into these problems 
emerges with the discovery of the Ai_yFe2-x^^2 fam- 
ily of FeSCs f8], where the considerable amount of 
Fe vacancies induce substantial changes in electronic and 
magnetic structures [9 . In particular, the parent com- 
pound K0.8Fe1.6Se2 exhibits an unusual insulating 2x2 
block checkerboard AF order (Fig. If, referred to X-type 
from now on) [TO-12 . The 20% Fe vacancies in it form a 
^/5 X ^/b order below T5'=578 K, on top of which the X- 
type spin order develops below T/v=559 K with a large 
ordered Fe magnetic moment 3.3 jUB- In a broad per- 
spective, K0.8Fe1.6Se2 has brought in an ideal bench- 
mark against magnetic theories for FeSCs: With fixed 
parameters good for vacancy-free chalcogenides, intro- 
ducing the ordered vacancies should transit the metallic 
E type to the insulating X type. In particular, under- 
standing the metal-insulation transition will yield insight 
into how electrons become correlated in FeSCs. 

So far, the weakly interacting itinerant-electron model 
based on the experimentally observed Fermi surface 
topology can reproduce the C type only |5l [T3HT5] . The 
Heisenberg spin model, as fit to neutron scattering data 
on CaFe2As2 [16 , FeTe [H], and K0.8Fe1.6Se2 [II], shows 
that the E and X types share similar model parameters. 



but the C and E types are surprisingly well separated 
in the model parameter space, with the leading exchange 
interaction being AF and ferromagnetic (FM), respec- 
tively. To recover the C — E proximity and their metal- 
licity, a model with coexisting itinerant electrons and lo- 
calized spins was proposed [18 . Interestingly, all these 
spin orders have been reproduced in first-principles band 
calculations [19-22 ; however, the microscopic origin of 
the magnetic softness and the metal- insulating transition 
remains to be elucidated. For example, the insulating 
nature of K0.8Fe1.6Se2 has been debated between Mott 
insulator [23] and magnetic semiconductor [21 scenar- 
ios. It is thus important and timely to carefully examine 
the electronic structure of vacancy-ordered K0.8Fe1.6Se2 
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FIG. 1. The in-plane patterns of the spin-up (blue balls) 
and spin-down (red balls) iron atoms in (a) collinear C-type, 
(b) bicollinear E'-type, and (c) block checkerboard X-type AF 
states. Their counterparts in the presence of x y/5 Fe va- 
cancy ordering are (d)-(f), respectively. Patterned horizontal 
and vertical solid lines represent dyz and dxz bonds, respec- 
tively; dashed lines denote the Fe vacancy broken bonds. 
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towards a unified picture. 

In this Letter, we present a combined first-principles 
and effective Hamiltonian analysis of the electronic struc- 
ture of vacancy-ordered K0.8Fe1.6Se2. Our first-principles 
results show (i) that its nonmagnetic Fermi surface does 
not display a nesting vector, and the bare spin suscepti- 
bility xo (q, = 0) is rather featureless and considerably 
weaker than the vacancy-free cases, rendering the Fermi 
surface instability an unlikely scenario, and (ii) that its 
magnetic electronic structure features two "gaps": a 
high-energy Mott gap and a low-energy Fe-Fe bonding 
gap, suggesting the picture of coexisting itinerant and 
localized electronic states [181 l24H27| . Consistently, we 
demonstrate that the insulating antiferromagnetism of 
K0.8Fe1.6Se2 can be unified with the metallic antifer- 
romagnetism of vacancy- free FeSCs in the spin-fermion 
model p^. These findings indicate that the blocking 
effects of Hund's rule coupling at low-energy scale and 
the resulting electronic correlation are crucial to the elec- 
tronic and magnetic structures of FeSCs [28} [29]. 

First-Principles Analysis [30 . — Since nonmagnetic 
K0.8Fe1.6Se2 is metallic [21], a first question that needs to 
be clarified is whether the magnetic and metal-insulator 
transitions in K0.8Fe1.6Se2 is driven by Fermi surface 
instability. To this end, the electronic structure of 
K0.8Fe1.6Se2 at the experimental crystal structure [31 
was calculated within local density approximation (LDA) 
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FIG. 2. LDA Fermi surface at kz=0 (a) and kz—'K (b) and 
band structure (c) in the one-Fe unit cell notation. High in- 
tensity means high spectral weight. Shadow bands appear as 
shifting the main bands by (tt, tt, tt) and Q = zb(37r/5, 7r/5, tt), 
=b(— 7r/5, 37r/5, tt). Color coded for dxz (red), dyz (blue), and 
the rest 3(i orbitals (green). Dashed lines in (a),(b) are the 
folded Brillouin zone boundary in the in-plane 8-Fe unit cell 
notation, (d) Bare spin susceptibility in the one-Fe unit cell 
notation, with the value at Q indicated by the arrow, com- 
pared with those of FeTe and a previous model [13] . 



of density functional theory (DFT) , implemented via full 
potential, all-electron, linearized augmented plane wave 
basis [32]. The resulting Fermi surface (more precisely 
the intensity of the one-particle propagator at Fermi 
level) unfolded [33 into the one-Fe unit cell notation is 
presented in Figs. 2(a) and 2(b) [34]. No Fermi surface 
nesting is visible at (37r/5,7r/5) and (— 7r/5, 37r/5), the 
characteristic wavevectors of the X-type spin order. 

To be quantitative, we calculated the bare spin sus- 
ceptibility xo(ct5^ = 0) and unfolded it into the one-Fe 
unit cell notation [30] . In comparison, we also calculated 
Xo(q7 0) for two vacancy- free cases having the Fermi sur- 
face nesting between hole pockets around (0, 0) and elec- 
tron pockets around (7r,0) and (0,7r): a one-Fe-unit-cell 
model [13] and FeTe in LDA. As shown in Fig. 2(d), the 
real part of Xo{^.->^) K0.8Fe1.6Se2 is rather featureless, 
in sharp contrast with the vacancy- free cases. Besides, its 
maximum value is substantially smaller, making it much 
harder for the interaction-assisted instability [13 to be 
effective. Hence, Fermi surface instability is unlikely a 
driving force for the magnetism in iron chalcogenides. 

Next, we elucidate the nature of the band gap opening 
in K0.8Fe1.6Se2. On the one hand, as shown in Fig. 2(c), 
removing one quarter of nearest neighbors in the Fe 
plane (Fig. If) reduces the bandwidth W of nonmagnetic 
K0.8Fe1.6Se2 from ~ 4 eV in vacancy- free FeSCs [34, 35 
to ~ 3 eV; thus, the larger U /W might favor the Mott 
metal-insulator transition, as illustrated in a two-orbital 
model without Hund's rule coupling [23^. On the other 
hand, the observation of a band gap A(/7 = 0) ~ 0.6 eV 
in generalized gradient approximation (GGA) of DFT, 
led to the conclusion that the X-type K0.8Fe1.6Se2 is a 
magnetic semiconductor [21]. Hence, it is critical to clar- 
ify the U dependence of A. 

To this end, we performed a series of GGA+/7 calcu- 
lations for a number of U on the Fe atoms, different spin 
orders, and Fe-Fe bond lengths, since a Mott gap gen- 
erally scales up with [/, independent of spin order, and 
is anticorrelated with the bond strength. As revealed in 
Fig. 3(a), the Fe partial density of states (DOS) shows 
that Fe 36^^ is always in the high-spin configuration (five 
spin- majority electrons and one spin- minority electron), 
manifesting the strong effect of Hund's rule coupling. It 
exhibits two "gaps" : One, a high-energy "Mott gap" be- 
tween the spin- majority subbands and spin-minority sub- 
bands, which scales with U and is insensitive to the spin 
order. Two, the low-energy real gap A resides at Fermi 
level (zero energy) within the spin-minority subbands. 

The fact that A resides in one spin channel suggests 
that a direct way to decern its nature is to manipulate the 
bond length within the FM 2x2 iron block (Fig. If), and 
that a key clue is the experimentally observed tetramer 
lattice distortion (TLD) [31 : The intra-block and inter- 
block Fe-Fe bond lengths are 2.691 and 2.916 A, respec- 
tively. We thus compared the X type in the realistic 
structure and in a hypothetic TLD-free structure where 
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FIG. 3. (a) Fe 3d partial DOS of the X-type AF state ob- 
tained from GGA+t/ calculations with [7 = and 3 eV. Up- 
per and lower panels are for spin-majority and spin-minority, 
respectively, (b) The band gap size as a function of U for 
X-type (solid squares) and FM (open circles) orders, as well 
as for the X-type without TLD (open triangles), (c) Total 
energy difference per Fe atom induced by vanishing TLD in 
the X type as a function of U. 
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FIG. 4. (a) Itinerant energy per Fe versus KS for the TLD 
parameter a = (see text). The FM state is the energy 
reference, (b) Itinerant energy per Fe versus a for KS = 0.8 
eV. (c) Gap size versus KS for a = 0. (d) Gap size versus a 
for KS = 0.8 eV. 



all the Fe-Fe bond lengths are equal to 2.757 A and the 
Fe-Se bond lengths remain unchanged j3Q]. 

We found that the Mott insulator scenario disfavors 
the realistic structure. The total energy difference be- 
tween these two structures as a function of U is plotted 
in Fig. 3c. The realistic structure turns out to be un- 
stable for > 3 eV. Consistently, as shown in Fig. 3b, 
the increase of A with U is relatively slow in the realistic 
structure for U < 3 eV (squares) but quick in the TLD- 
free structure for /7 > 3 eV (triangles). Moreover, in the 
realistic regime {U < 3 eV), A in the TLD-free structure 
is noticeably smaller than in the realistic case, indicating 
a positive correlation between A and intra-block Fe-Fe 
bond strength. Furthermore, A is especially sensitive to 
the magnetic structure for /7 < 3 eV, as manifested by its 
vanishing values in the FM case (Fig. 3b). Hence, the re- 
alistic A results essentially from the bonding-antibonding 
splitting within the 2x2 FM iron block. 

Effective Hamiltonian. — ^The above first-principles re- 
sults suggest the picture of low-energy itinerant electrons 
spin-polarized by Hund's rule coupling to more localized 
electronic states. A minimum model is the spin-fermion 
model where the Fe dxz and dyz orbitals were treated to 
host itinerant electrons and the rest Fe 3d orbitals were 
treated as an effective localized spin [181 ESI ES] • 

where C^^^ denotes the annihilation operator of an itiner- 
ant electron with spin ja =t or ^ in the 7 = d^z or dy^ or- 
bital on site i. tj^ 's are the electron hopping parameters. 



d^n' is the Pauli matrix and Si is the localized spin whose 
magnitude is S. K is the effective Hund's rule coupling. 
Jij is the AF superexchange couplings between the lo- 
calized spins; in particular, J and J' are respectively the 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
ones. The filling of the itinerant electrons is on average 
three (one hole) per Fe site, corresponding to the high- 
spin configuration of Fe 3d^ [35^ . 

This model was proposed [18] to unify the metallic 
C-type and £'-type AF orders in vacancy-free FeSCs, as- 
suming that KS is the leading material-dependent pa- 
rameter, controlled by the anion height from the iron 
plane; KS was set to be 0.8 eV for FeTe. To show the 
unifying capability of the model, below the same set of 
parameters is used for K0.8Fe1.6Se2. The effect of TLD 
in K0.8Fe1.6Se2 is included via multiplying the intrablock 
(and interblock) hopping parameters by 1+a (and 1 — a), 
corresponding to the shrinking (and elongation) of the 
Fe-Fe bonds. Comparing the hopping parameters ob- 
tained from our first-principles Wannier function analy- 
sis [3F, we found a 0.2. The magnetic landscaping is 
studied by comparing a variety of static spin orders in 
the presence of ordered vacancies, such as the FM state 
and the AF states of C-type (Fig. Id), E'-type (Fig. le), 
X-type (Fig. If), and G-type (i.e., the Neel state where 
all NN spins are antiparallel) , with the localized spins 
treated as Ising spins. 

Figs. 4a and 4b show how the X type is stabilized. The 
localized-spin part of the model favors the C type due to 
the comparable J and J' superexchange processes: Their 
contributions to the total energy per iron is —l.bS'^J^ for 
the C type and 0.55^(7— J') for the X type, for example. 
It is thus the itinerant energy, the energy involving the 
itinerant electrons, that favors the X type around KS = 
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0.8 eV (Fig. 4a). The itinerant energy is further lowered 
by TLD a remarkably (Fig. 4b) enough to make the X 
type more stable than the C type for J' S'^ < 30 meV. 
These results agree well with the earlier ah initio study 
of the TLD effect [21^. Furthermore, in the vicinity of 
the X and C types, the FM state is considerably higher 
in total energy than all these AF states by more than 
50 meV per Fe atom (not shown), in agreement with 
the first-principles total energy calculations [211 El], a 
feature current spin-only model analysis failed to capture 
[TTJ[22]. These results imply that with moderate KS^ this 
system warrants strong and overall AF spin fluctuations 
in the Fe plane, providing a necessary environment for 
singlet superconductive electron pairing. 

The insulating property of the X type is also repro- 
duced. The gap size as a function of KS and a is shown in 
Figs. 4c and 4d, respectively. The nature of the band gap 
can be inferred from the limit of large KS^ the so-called 
double-exchange limit where the electrons cannot hop to 
the sites with opposite spin orientation due to the en- 
ergy barrier as high as KS. In this limit, the inter-block 
hoppings occur at distances longer than NNN (Fig. If). 
Thus, with the NN and NNN hopping parameters, any 
itinerant electron is localized within a 2 x 2 block, lead- 
ing to four discrete energy levels. Even for KS = 0.8 eV, 
the inter-block hoppings are still strongly suppressed and 
cannot help develop overlapped bands from those energy 
levels, as predicted in Ref. jl^ (FeTe is metallic because 
of unsuppressed NNN hoppings; see Fig. lb). Therefore, 
the band gap originates from bonding-antibonding split- 
ting in the spin-minority channel and its size is positively 
correlated with the intra-block Fe-Fe bonding strength 
represented by a, indeed. 

The above recognition that optimization of the Fe-Fe 
FM bonds is an important factor leads to the following 
intuitive insight into the vacancy- induced E-X transi- 
tion. As shown in Figs, lb and If, both spin orders 
share a similar bond pattern: Each Fe atom is linked 
to one dxz FM bond and one dyz FM bond. We notice 
that another realization of this pattern is the vacancy-free 
X type (a metallic 2x2 block checkerboard AF order) 
shown in Fig. Ic. Thus, the vacancy-free E and X types 
are likely to have similar itinerant energy. Besides, their 
localized-spin part contributes exactly the same energy 
in the classic spin approximation. This implies that the 
vacancy-free X and E types could be very close in en- 
ergy, in agreement with neutron scatter measurement on 
FeTe [36]. Indeed, we confirmed that they differ by only 
7 meV/Fe in ah initio calculations and 3 meV/Fe in Eq. 
(1) for FeTe. Then, the added Fe vacancies break the 
bonding pattern in the E type (Fig. le) but retain it in 
the rearranged X type (Fig. If). This means that in iron 
chalcogenides, the 2x2 block checkerboard AF order is 
already highly competitive in the absence of Fe vacancies, 
and emerges as the ground state upon introduction of Fe 
vacancies. This can also serve as the base to understand 



the insulating 2x2 block AF order in BaFe2Se3 [37] , 

It is noteworthy that for a system with degenerate or- 
bitals, orbital ordering may occur. On-site orbital order- 
ing was argued to drive the C type in LaOFe As |^ , the 
E type in FeTe [38^ and the X type in Ko.sFei.eSes [39]. 
The aforementioned peculiar bond pattern implies that 
the on-site orbital polarization P (the difference in the 
occupation numbers of the Fe dxz and dyz orbitals) is 
weak in K0.8Fe1.6Se2. Our first-principles Wannier func- 
tion analysis verified that P ~ in FeTe and P 0.06 in 
K0.8Fe1.6Se2, much smaller than P = 0.17 in LaOFe As 
[35 . The nonvanishing P in K0.8Fe1.6Se2 is caused by the 
vacancy-imposed symmetry breaking, whose effect is sur- 
prisingly weak, confirming that the Fe-Fe bonding effect 
is dominant. Interestingly, while on-site orbital ordering 
is found weak (confirmed in Ref. [40]), the bond pattern 
may be regarded as a bond orbital order. 

The importance of the Hund's rule coupling K in 
governing electronic correlation in the metallic state of 
vacancy- free FeSCs has also been recently demonstrated 
in LDA and dynamical mean-field theory, giving rise to 
a new term "Hund's metal" [28l[29]. Likewise, insulating 
K0.8Fe1.6Se2 may be regarded as a realization of novel 
"Hund's insulator" in the sense that the band gap opens 
primarily to gain Hund's rule coupling energy rather than 
Hubbard repulsive energy (i.e., due to the blocking effect 
of K rather than U). A measure to distinguish Hund's 
insulator and Mott insulator is checking whether the gap 
size is sensitive to the magnetic structure rather than 
U or positively correlated with the FM bond strength. 
Thus, Hund's metal-insulator transition is expected to be 
much more sensitive to structural changes. This suggests 
that tuning the Ai-yFe2-xSe2 materials through Hund's 
metal-insulator transition (e.g, by pressure) be an effec- 
tive route to optimize their superconductive properties. 

In summary, we have shown from first principles that 
the insulating X-type antiferromagnetism in vacancy- 
ordered K0.8Fe1.6Se2 is incompatible with a simple 
Fermi-surface nesting or Mott insulator scenario. We 
demonstrate that it can be unified with the metallic 
C-type and £^-type antiferromagnetism of vacancy-free 
FeSCs in the spin-fermion model, where the competi- 
tion between double-exchange ferromagnetism associated 
with Hund's rule coupling and superexchange antiferro- 
magnetism is tuned by the ordered vacancies into a novel 
Hund's metal-insulator transition. These findings indi- 
cate that the blocking effects of Hund's rule coupling and 
the resulting electron correlation are crucial to the elec- 
tronic and magnetic structures of FeSCs, and are likely 
at the heart of their high-Tc mechanism. 
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